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The distance and internal composition of the neutron star 
in EXO 0748-676 with XMM-Newton 
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ABSTRACT 

Recently, the neutron star X-ray binary EXO 0748-676 underwent a transition to 
quiescence. We analyzed an XMM-Newton observation of this source in quiescence, 
where we fitted the spectrum with two different neutron-star atmosphere models. From 
the fits we constrained the allowed parameter space in the mass-radius diagram for 
this source for an assumed range of distances to the system. Comparing the results 
with different neutron-star equations of state, we constrained the distance to EXO 
0748-676. We found that the EOS model 'SQM1' is rejected by the atmosphere model 
fits for the known distance, and the 'AP3' and 'MSI' is fully consistent with the known 
distance. 
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1 INTRODUCTION 

The low-mass X-ray binary (LMXB) EXO 0748-676 was 
discovered as a transient source with the European X- 
ray O bservatory Satellite (EXOSAT) in 1985 (|Parmar et all 
1986). The source exhibits simultaneous X-ray and opti- 
cal eclipses fr om which the orbital period of P = 3.82 hr 
was deduced l|Crampton et al.lll986t ). EXO 0748-676 also 
exhib ited irregular X-ray dipping activity (IParmar et al.l 
Il986h . and type-I X-ray bursts (JGottwald et al.l 1 19861) . 
Burst oscillati ons in EXO 0748— 676 were first reported 
by I Villarreal &: Strohmaveri |2004) at 45 Hz in the average 
Fourier Power Spectrum of 38 type-I X-ray bursts; the 45-Hz 
signal was then int erpreted as the spin fre quency of the neu- 
tron star. Recently. iGallowav et al.l (|2009l ) detected millisec- 
ond oscillations in the rising phase of two type-I X-ray bursts 
in EXO 0748-676 at a frequency of 552 Hz. They concluded 
that the spin frequency of EXO 07 48-676 is close to 522 Hz 
rather than 45 Hz as suggested by I Villarreal fc Strohmaveri 
(2004). The 45 Hz oscillation may arise in the boundary 
layer between the disk and the neutr on star (Balman 2009) 



or it could be a statistical fluctuati on (|Gallowav et al.ll2009l ) 



ICottam. Paerels fc Mendez! (|2002l ) reported a measurement 
of the gravitational redshift from iron and oxygen X-ray 
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absorption lines arising from the atmosphere of the neu- 
tron star in EXO 0748—676 during type-I X-ray bursts, 
but subsequent obse rvations failed to confirm these features 
llCottam et aU 120081. Based on the gravitational redshift, 
IQzell (|2006 ) suggested that the mass, radius and distance of 
EXO 0748-676 are 2.10±0.28 M , 13.8±1.8 km and 9.2±1.0 
kpc, respectively, which would rule out many neutron-star 
equations of state. 

Measuring the distance to LMXBs is difficult, except 
for the case of sources in globular clusters. A way to get 
the distance is using type-I X-ray bursts. The peak flux 
for some very bright bursts can reach the Eddington lu- 
minosity at the surface of the neutron star. From a strong 
X-ray burst, I Wolff et at] (|2005l l derived a distance to EXO 
0748—676 of 7.7 kpc for a helium-dominated burst photo- 
sphere, and 5.9 kpc for a hydro gen-dominated burst photo- 
sphere. IGallowav et al.l (|2008al l analyzed several type-I X- 
ray bursts from EXO 0748—676 and estimated a distance 
of 7. 4 kpc, different from the value of 9.2 kpc reported by 
IQzell (|2006l ). Taking into account the touc hdown flux and 
high-in clination in EXO 0748—676 , recentlv lGallowav et al.l 
(2008b) gave a distance of 7.1 ± 1.2 kpc. 

Another way to get the distance to an LMXB is through 
observations of quiescent X-ray emission from the neutron- 
star surface. During the quiescent state, X-ray emission orig- 
inates from the atmosphere of the neutron star. By fitting 
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the X-ray spectrum of the neutron-star system with hydro- 
gen atmosphere models, one can estimate the mass, radius 
and distance of the neutron star. Recently, the neutron-star 
X-ray trans ient EXO 0748—676 underwent a transit ion into 
quiescence (jDegenaar et alj|2009l ; iBassa et al.ll2009T ). 

In this paper, we report on the distance to EXO 
0748-676 that we constrained from XMM-Newton data. We 
use two different neutron-star atmophere models to fit the 
X-ray spectrum, and compare the results of the spectral fit- 
ting with different neutron-star equation of state (EOS). In 
the next section, we describe the observation and data anal- 
ysis. We show the fitting results in §3l and we discuss our 
findings in §3] 



2 OBSERVATIONS AND DATA ANALYSIS 

EXO 0748—676 was observed with the European Photon 
Imaging Camera (EPIC PN and MOS) on board the XMM- 
Newton on 2008 November 6 at 08:30:03 UTC (obsID 
0560180701). The PN and the two MOS cameras were oper- 
ated in Full- Window mode. We reduced the XMM-Newton 
Observation Data Files (ODF) using version 8.0.0 of the sci- 
ence analysis software (SAS). We used the epproc and em- 
PROC tasks to extract the event files for the PN and the two 
MOS cameras, respectively. Source light curves and spectra 
were extracted in the 0.2 — 12.0 keV band using a circular 
extraction region with a radius of 30 arcsec centered on the 
position of the source. Background light curves and spec- 
tra were extracted from a circular source-free region of 35 
arcsec source- free on the same CCD. We applied standard 
filtering and examined the light curves for background flares. 
No flares were present and we used the whole exposure for 
our analysis. The exposure time for the PN camera was 24.2 
ks, and for each MOS camera was 29.03 ks. The source count 
rate was 0.496 ± 0.005 cts/s for PN, and 0.135 ± 0.002 cts/s 
and 0.127 ± 0.002 cts/s for MOS1 and MOS2, respectively. 
We checked the filtered event files for photon pile-up by run- 
ning the task epatplot. No pile-up was apparent in the PN, 
MOS1 and MOS2 data. The photon redistribution matrices 
and ancillary files for the source spectra were created us- 
ing the SAS tools rmfgen and arfgen, respectively. We 
rebinned the source spectra using the tool PHARBr>Q such 
that the number of bins per resolution element of the PN 
and MOS spectra was 3 and the minimum number of counts 
per channel was 20. 

We fitted the PN and MOS spectra simultaneously 
in th e 0.5-10.0 keV range with XSPEC 12.50 (JAmaudl 
|l996J), using either of t wo neutron-star hyd rogen-atmosphere 
models: NSAGRAV i|Zavlin et al.l 1 1996ft and NSATMOS 
l|Heinkeet all 120061 ) . The NSAGRAV model provides the 
spectra emitted from a nonmagnetic hydrogen atmosphere 
of a neutron star with surface gravitational acceleration, g, 
ranging from 10 13 to 10 15 cm s~ 2 . This model uses the mass 
(Mns) and radius (-Rns) of the neutron star and the unred- 
shifted effective temperature of the surface of the star (kT e g) 
as parameters. The normalization of the model is defined as 
1/.D 2 , where D is the distance to the source in pc. The 
second model that we used, NSATMOS, includes a range 
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Figure 1. XMM-Newton PN (black), MOS1 (red) and MOS2 
(green) spectrum of EXO 0748-676 in the 0.5- 10.0 keV en- 
ergy band. The spectrum was fitted with a neutron-star hydro- 
gen atmosphere model (NSATMOS) and a power-law model with 
r fixed to 1. The lower panel shows the residuals to the best-fit 
model. 



of surface gravities and effective temperatures, and incor- 
porates thermal electron conduction and self-irradiation by 
photons from the compact object. This model assumes neg- 
ligible magnetic fields (less than 10 9 G) and a pure hydrogen 
atmosphere. NSATMOS parameters are Mns, -Rns, logT e g 
(the same as for NSAGRAV), distance in kpc, and a sepa- 
rate normalization K, which corresponds to the fraction of 
the neutron-star surface that is emitting. We fixed K to be 
1 in all our fits with NSATMOS. 

We included the effect of interstellar ab- 
sorption using PHABS assuming cross-sections of 
I Balucinska- Church fc McCammonl Jl992j) and solar abun- 
dances from lAnders fc Grevessd i|1989l) . and we let JVh, 
column density along the line of sight free to vary during the 
fitting. In order to account for differences in effective area 
between the different cameras, we introduced a multiplica- 
tive factor in our model. First, this factor was fixed to unity 
for PN and free for MOS1 and MOS2. Then, we set the 
scaling factor to unity for MOS1 and MOS2, respectively, 
and set the factor free for the other cameras. We found 
that, fixing the scaling factor for different cameras gives 
similar best-fit results. Therefore in the rest of the paper 
we fixed the factor to be 1 for PN and free to vary for 
the other cameras. None of the atmosphere models alone 
fitted the spectrum above ~ 2 — 3 keV properly. Adding 
a power-law component improved the fits significantly, 
however, all parameters were less constrained than when 
fitting the data with the neutron-star atmosphere model 
only. We first fixed the power-law index to 0.5, 1.0 and 
1.5 to get better constraints on the paranreters of the 
neutron-star atmosphere nrodel (Degenaar et al. 2009). 
Further, we initially fixed the distance to the NS at 7.1 
kpc, which is the value inferred from the touchdown flux of 
iGallowav et"afl (|2008bl ) . 
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Table 1. Best-fit parameters of neutron-star atmosphere models fit to the XMM-Ncwton data of EXO 0748—676 



model 


N H 
(10 20 em- 2 ) 


J off 

(eV) 


M NS 
(M©) 


Kns 
(km) 


r 


10 — 13 ergs cm -2 s _1 


Fx 
10 -12 ergs cm -2 s~ 1 


X 2 /d.o.f. 


NSAGRAV 


5.6 ±1.8 


U3^ 4 


1.55 ±0.18 


15.2 ±1.8 


0.5 


1.15 ±0.21 


1.18 ±0.15 


0.986/219 


NSATMOS 


5.4 ±1.5 


113 ±4 


1.29 ±0.20 


16.11S1 


0.5 


1.17 ±0.20 


1.23 ±0.16 


0.985/219 


NSAGRAV 


6 2+ 13 


114±f 


1.62 ±0.11 


15.8i°f 


1.0 


1.10 ±0.15 


1.14 ±0.13 


0.977/219 


NSATMOS 


6.1 ±1.5 


114 ±4 


1.55 ±0.12 


i6.ot; : ; 


1.0 


1.11 ±0.15 


1.13 ±0.06 


0.977/219 


NSAGRAV 


6.7 ±1.5 


110 ±8 


1.71 ±0.30 


16.5 ±0.5 


1.5 


1.00 ±0.19 


1.01 ±0.15 


0.987/219 


NSATMOS 


6.7 ±1.4 


110 ±5 


1.77 ±0.45 


16 -«±T.B 


1.5 


1.03 ±0.22 


1.03 ±0.10 


0.985/219 



Note. — Njj is the equivalent hydrogen column density, T^ the effective temperature of the neutron-star surface as seen at infinity, 
MjMg and -Rims are the mass and radius of the neutron star, respectively. F pow is the unabsorbed flux of the power-law component in 
the 0.5—10 keV energy band, and Fx is the total unabsorbed X-ray flux in the same energy band. The last column gives the reduced 
X 2 for 219 degrees of freedom. The quoted errors represent the 90% confidence levels. 



3 RESULTS 

3.1 Results from the spectral fits 

Figure [Q shows the XMM-Newton spectra of EXO 0748-676 
fitted with the model "phabs (NSATMOS + powerlaw) ". 
The power-law index is fixed at 1.0. The best fit of this model 
gives -/Vh = 6.1 ±1.5 xlO 20 cm -2 , neutron-star mass Mns = 
1.55±0.12M Q , neutron-star radius _Rns — 16.0+5'g km, and 
effective temperature logoff = 6.20 ± 0.02 (in K). Accord- 
ing to th e same formula Tig = T e gi/l - (2GM N s)/(-Bnsc 2 ) 
used by iDegenaar et al.l l|2009l ). we converted T e g to the 
effective temperature as seen by an observer at infinity, 
T£§ = 114 ±4 eV. In the formula, G is the gravitational con- 
stant and c is the speed of light. The model predicts 0.5—10 
keV an unabsorbed X-ray flux F x = 1.13 ± 0.06 x 10~ 12 
ergs cm -2 s _1 . The flux of the power-law component in the 
same energy band is F pow = 1.11 ± 0.15 x 10 -13 ergs cm -2 
s _1 , which corresponds to ~ 10% of the total unabsorbed 
flux. The reduced \ 2 is 0.977 for 219 degrees of freedom. 
The best-fit results of the models NSAGRAV and NSAT- 
MOS for the three different power-law index are given in 
Table-[T] Errors are given at the 90% confidence level for one 
fit parameter. 

We note from Table \T\ that both atmosphere models, 
regardless of the value of F, yield a good fit with similar 
X 2 ■ In the rest of the analysis, we used a power-law index 
fixed to 1. Further, A?h and T e g are well constrained and 
are consistent for the different fits. Both NSAGRAV and 
NSATMOS models also give consistent results on Mns and 
.Rns- The NSATMOS model is more accurate in constraining 
T cff than the NSAGRAV model. 



3.2 Equation of state 

Fitting the quiescence XMM-Newton spectrum of EXO 
0748—676 with two different atmosphere models and com- 
paring the results allows us to test the reliability and ac- 
curacy of both models. From the fits we get a mass and 
radius of the neutron star at a specified distance, and then 



by comparing the inferred mass and radius with the differ- 
ent neutron-star EOS we can give upper limits to the source 
distance for the different EOS. 

We used the steppar command in xspec to vary the 
mass, radius and distance parameters simultaneously, allow- 
ing other parameters to be free to find the best fit at each 
step. For the mass we go from 0.5 to 2.5 Mq with steps of 0.1 
Mq, and for the distance we go from 5 to 10 kpc with steps 
of 0.25 kpc. The minimum and maximum radius allowed 
with these models are 5.0 km and 25.0 km, respectively. In 
Fig [2] we show the contour plots obtained from the STEP- 
PAR procedure for the NSATMOS model. Each plot is for 
a different distance, ranging from 5 to 10 kpc. The contour 
lines (red) are for the confidence levels of 90% (solid) and 
99% (dashed). Further, in Fig [2] we give different neutron - 
star EOS (black) taken from lLattimer fc Prakashl (|2007l ). 
We did the same analysis for the NSAGRAV model as well. 
Both two mod els give consistent resu lt, in accordance with 
the findings of I Webb fc Barretl (|2007l ). 

Using the optical data from the Very Large Tele- 
scope (VLT), moderate-resolution spectrosco py of the op- 
tical counterpart and Doppler tomography, iMurioz et al.l 
(2009) provided the first dynamical constraints on the stel- 
lar mass of LMXB EXO 0748-676 . The mass range of the 
neutron star that they de r ived i s IMq ^ Mns ^ 2.4M©. 
Subsequently, Bas sa et al.l (|2009l ) analyzed optical spectra 
of EXO 0748—676 when the source was in the quiescent 
state, and they gave a lower limit to the neutron-star mass of 
M ns ^ 1-27 Mq . As up per limit we used the value reported 
by iMunoz et al.l l|2009l ). but since at the time of their ob- 
servation the so urce was still in out burst, we used the lower 
limit reported in lBassa et all (2009). In Figure[2]we also give 
the lower (pink/dotted) and upper (green/dashed) limits to 
the neutron-star mass. 

In order to test the EOS and identify the upper limit to 
the source distance, we assume three different EOS models: 
normal nucleonic matter (AP3), boson condensates matter 
(MSI) and strange quark matter (SQM1). By varying the 
source distance from 5 to 10 kpc, the contour lines for the 
fitted model move on the NS mass-radius diagram. We can 



4 G. Zhang et al. 





2.5 




^^^p: 


JPA1 


msAmso 


Munos 09 


90% 
99% 







■^-FNgN:-"^ 


v I 












2 




SQM3Jteu\p 


A' \ 

vMS1 \ 












1.5 


Z-— SQMJ^ 

""" gsTVT 


^PALe h~ 


3M3k 






Bassa 09 




" 




1 
0.5 




il 






£r-- 



10 12 14 16 18 20 

Radius (km) 

(a) 5.0 kpc 



2.5 


^F" ""N. 90 % — ' 

^^^ MS? \MSO 
^R^-v^lPA i \ 09 


2 

1 1.5 

1 

0.5 


■^^ !-"--IC 
Bj^^ AP4V\ 

S Q, M1 PAL6 % 
GS1 ] " --^/ 


\ \ NAPL1 
3UTt4-sMS1 \ 


Bassa 09 



10 12 14 16 18 

Radius (km) 

(b) 6.0 kpc 




10 12 14 16 

Radius (km) 

(c) 7.0 kpc 



8 20 




10 12 14 16 18 20 

Radius (km) 

(d) 8.0 kpc 




^^^^^^^^^^■P" "~^\ 90 % — 

^^^ MS2 \ MS " : ' »% 

^■fff^Jill-'AI NK * : \ Munos 09 j 


^^ SQM3J^ 
SQM1 pA| _ 6 / 

gsij — ^->/ 


\ V"^PL1 

^3r\J >, \ 


Bassa 09 '., 


-^/ 




l "1 




fe ;/ : 



Radius (km) 

(e) 9.0 kpc 



Radius (km) 

(f) 10.0 kpc 



Figure 2. Contour plots showing the results of modeling the neutron-star in EXO 0748—676 with the XSPEC model NSATMOS and 
power-law. The power-law index is fixed to 1. The plots show two confidence levels in the mass-radius diagram obtained from our fit; 
the contour lines (red) are for the confidence levels of 90% and 99%, respectively. The pink line "A" is for the lower limit of Mjms given 
by Bassa et al.(2009), and the green line "B" is for the upper limit given by Muiioz et al. (2009). 



estimate the probability of the distance for each EOS when 
the contour lines pass through the EOS curves. Note that 
as the distance increases (see Figure [2]), the satisfied area of 
the model moves from bottom left to top right in the plot. 
The results using NSAGRAV are similar to those shown in 
Figure [2] For a certain distance we found that not all the 
EOSs are consistent with the two neutron-star atmosphere 
models that we used. 

If the neutron star in EXO 0748-676 follows the EOS 
model 'AP3', the probability that the source has a distance 
of 10.0 kpc is 1 x 10~ 4 and 1 x 10~ 6 for NSAGRAV and 
NSATMOS, respectively. If we want to get a probability 
for the distance larger than 1 x 10 -2 (99% confidence), the 
distance for NSAGRAV and NSATMOS should be smaller 



than 8.9 kpc and 8.5 kpc, respectively. The distance at 90% 
confidence for NSAGRAV and NSATMOS is less than 8.3 
kpc and 8.2 kpc, respectively. Both models are consistent 
with t he distan ce of 7 .1 kpc given by type-I X-ray bursts 
l|Gallowav et al.ll2008bh . 

For the EOS model 'MSI', the probabilities that EXO 
0748-676 is at a distance of 10 kpc is 10~ 5 and 10~ 6 for 
NSAGRAV and NSATMOS, respectively. For both models, 
respectively, the distance at 99% confidence level is less than 
7.3 kpc and 7.1 kpc, and the distance at 90% confidence 
level is less than 6.9 kpc and 6.8 kpc. Both neutron-star 
atmosphere models with the 'MSI' model have an upper 
limit for the distance smaller than 7.1 kpc. 

For a EOS model 'SQM1', the distance at 99% confi- 
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Table 2. Upper limits on the distance to EXO 0748-676 for 
different EOS models. 



EOS 


AP3 


AP3 


MSI 


MSI 


SQM1 


SQM1 


confidence 


90% 


99% 


90% 


99% 


90% 


99% 


NSAGRAV 


< 8.3 


< 8.9 


< 6.9 


< 7.3 


< 5.0 


< 5.2 


NSATMOS 


< 8.2 


< 8.5 


< 6.8 


< 7.1 


< 5.0 


< 5.2 



Note. -The 90% and 99% confidence levels upper limit for the 
two NS atmosphere models NSAGRAV and NSATMOS for 
the EOS models: 'AP3', 'MSI' and 'SQM1'. The distance is 
in kpc. 

dence level is less than 5.2 kpc, and the distance at 90% 
confidence level is less than 5.0 kpc for both atmosphere 
models. The upper limits on the distance to EXO 0748—676 
for different EOS are shown in Tabled The 'SQMT model 
is rejected at a 99% confidence level for this neutron star, 
unless the source is closer than 5.2 kpc. 



the NSAGRAV model are slightly higher than those for the 
NSATMOS model. The EOS model 'MSI' can be just satis- 
fied at a distance of 7.1 kpc. If we assume that the neutron 
star in EXO 0748—676 is a normal neutron star, following 
the EOS 'AP3', the source should be closer than 8.9 kpc for 
the NSAGRAV model, or 8.5 kpc for the NSATMOS model. 
Both the 'MSI' and 'AP3' EOS ar e fully consistent with 
the measured di stance of 7.1 kpc (|Gallowav et al.l l2008bl ; 
IWolff et atj|2005f ). For larger distances more EOS are ruled 
out. The EOS 'SQMT is rejected by the atmosphere model 
fits for a distance of 7.1 kpc measured from the X-ray bursts 
iGallowav et al.l (|2008bl ). We note, however, the neutron-star 
atmosphere models may not appropriate for 'bare' quark 
matter stars, but only for those normal quark star where a 
crust is present. 



ACKNOWLEDGMENTS 

This work is based on the observations obtained from XMM- 
Newton. PGJ acknowledges support from a VIDI grant from 
the Netherlands Organisation for Scientific Research. 



4 DISCUSSION 

We analyzed an XMM-Newton observation of the neutron 
star EXO 0748—676 in the quiescent state. The unab- 
sorbed X-ray flux in the 0.5—10.0 keV energy band was 
~ l.lx 10~ 12 ergs cm -2 s _1 . We found that the non- 
thermal (power-law) component only contributes ~ 10 ± 2% 
of the 0.5—10 keV X -ray flux, which is lower than what 
iDegenaar et al.l (J2009J) found from Chandra data (F pow was 
- 16-17% of the 0.5-10 keV X-ray flux from the fit with 
r = 1) about a month earlier than our observation. The total 
unabsorbed flux (0.5—10.0 keV) decreased from 1.3 x 10~ 12 
ergs cm -2 s _1 in the Chandra observation to 1.1 x 10~ 12 
ergs cm" 2 s _1 in our observation, whereas Nh changed 
from ~ 1.2 xlO 21 cm -2 to ~ 0.6 xlO 21 cm -2 . The effec- 
tive temperature, however, did not show large variations in 
one month time. According to the above comparisons, the 
reduction of the total flux is due to a lower contribution of 
the power- law component. 

Because the X-ray spectrum in the quiescent state is 
dominated by thermal emission originating from the NS sur- 
face, our data allow us to constrain the mass and radius 
of the neutron star. From the two different NS atmosphere 
models (NSAGRAV and NSATMOS) that we used to fit 
the X-ray spectrum, we found that both models show simi- 
lar results and set good constraints on the neutron-star ra- 
dius. Even takin g into account the Mns lower limit ( from 
iBassa et afll2009h . upper limit ffrom lMunoz et alj[200gh and 
our best fit A^ 2 contour, we still have a large area on the 
mass-radius diagram, and many EOSs are still possible (see 
Figure [2J) • In order to constrain the allowed space of mass 
and radius at a specified distance, we choose three typical 
neutron-star EOS, 'AP3', 'MSI' and 'SQMT. We found that 
the smaller the distance to the NS the more EOSs are con- 
sistent with the data. 

For any specific EOS, as the upper limit of the dis- 
tance we took the value of the distance where the 99% con- 
fidence contour just intersects the curve of that EOS. We 
found that the upper limits on the distance as derived from 
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